International Journal of Advanced Engineering Research 


and Science (IJAERS) 


Peer-Reviewed Journal 


ISSN: 2349-6495(P) | 2456-1908(O) 
Vol-9, Issue-4; Apr, 2022 
Journal Home Page Available: https://ijaers.com/ 

Article DOI: https://dx.doi.org/10.22161/ijaers.94.40 


International Journal of Advanced 
Engineering Research and Science 





Analysis nonlinear vibrations of the three - phase 
composite shallow cylindrical shell 


Huyen TTH, Luot QT 


Department of Mechanics, Le Quy Don Technical University, Hanoi City 100000, Vietnam 


Email: huonghuyen.hvktgs @ gmail.com 


Received: 21 Mar 2022, 

Received in revised form: 20 Apr 2022, 
Accepted: 25 Apr 2022, 

Available online: 30 Apr 2022 


©2022 The Author(s). Published by AI 
Publication. This is an open access article 
under the CC BY license 


Abstract— This paper presents the results 


the study nonlinear 


vibrations of the three-phase composite shallow cylindrical shell. The 
differential equation describing the nonlinear vibrations of the shell is 
solved by the Newmark direct integration method combined with the 
Newton Raphson iterative method. Numerical simulation by finite 
element method is used to calculate the vibration of the structure. The 
results of the digital survey allow to make quantitative comments, 
technical recommendations, help the designers and users to orient 


(https://creativecommons.org/licenses/by/4.0/). | effective applications in technical fields. 


Keywords— Nonlinear vibrations; Three- 
phase composites. 





I. INTRODUCTION 


Three-phase composite is a material consisting of a 
matrix phase, a fiber phase and a particle phase, which has 
been studied by Vanin G.A and Duc N.D. (1996a, 
1996b)[7,8]. In the 1997 publications, Vanin G.A. and Duc 
N.D determined the elastic modulus for three-phase 
composite materials 3Dm (Vanin G.A. and Duc N.D., 
1997) and 4Dm (Duc N.D., 1997a). An overview of three- 
phase composite materials was also found in the study of 
Minh D.K. (2011). Recently, Duc N.D. et al. (2011) also 
studied the nonlinear stability of three-phase polymer 
composite panels under thermal and mechanical load 
conditions (Duc N.D. et al., 2013; 2014). In this report, the 
authors studied the nonlinear dynamic response of the 
multilayer three-phase composite shallow cylindrical shell. 
The formulas are based on classical shell theory, taking 
into account geometric nonlinearities. 


I. MATHEMATICAL MODELS AND 
ASSUMPTIONS 


The cylindrical shell is considered as the shell with the 


ratio T < 5 , Where fo - the curvature of the shell, lmin = 


min 
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min(a, L), a, L- the equal projection dimensions of the 
shell (Fig.2). According to the finite element method, the 
cylindrical shell can be discretized by flat elements, 
whereby the shell is a finite combination of 9 node-pointed 
flat elements, called the flat shell element, where each “flat 
shell element” can be seen as a combination of two types 
of elements: a 9 node flat element, each with 2 degrees of 
freedom (ui, vi) and a 9 node flat shell element subjected to 
combined bending - torsion, each node has 4 degrees of 
freedom (Wi, Oxi, Oyi, Ozi), as shown in Fig.2. In the report, 
the authors use the thick-shell theory, which satisfies the 
Reissner-Mindlin theory. 














a, Shell model, b, Three-phase polymer composites model 


Fig. 1(a,b): The model of the problem 
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b, Flat element in tension (compression)and flat shell 
element in combined bending - torsion 


Fig. 2 (a,b): Finite Element Model 


HI. DETERMINATION OF ELASTIC 
COEFFICIENTS FOR THREE-PHASE 
COMPOSITE 


For three-phase composite materials (polymer matrix, 
fiber and particle). According to [2], the elastic coefficient 
of three-phase composites Vanin G.A. determined in 2 
Steps: 


The first step: Considering a two-phase composite 
consisting of the initial matrix phase and particles, such a 
composite is considered to be homogeneous, isotropic, and 
has two elastic coefficients. The determination of the 
elastic coefficients for composites filled with spherical 
particles is determined, taking into account the interaction 
between the particles and the matrix. The elastic 
coefficients of the grain-reinforced composite are now 
called hypothetical composites. 


Second step: determine the elastic coefficients of the 
composite between the assumed foundation and the 
reinforcing fibers. 


Assuming the components of the composite (matrix, 
fibers, particles) are all homogeneous, isotropic, then we 
will denote: En, Gn, Vn, Yn; Es, Gs, Vs, Ws; En, Gn, Vh, Wh are 
denoted by the modulus of elasticity, modulus of elasticity 
of shear, modulus of volume deformation, Poisson's 
coefficient, and composition ratio (by volume) of the 
matrix and particles, respectively. From here on, the 
quantities related to the matrix will have the n-index; 
relative to the particle is the h-index. According to [3], the 
elastic modulus of the assumed composite as follows: 
9KG .ṣ K-20 (1) 


E=—=; ———= 
3K +G 6K -2G 
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Where: 


— l4 A, G,LG K, yt (2) 
"1+y,(8—-10v,)H ` 


"1-4y, G.LGK,)" 


with 
K, -K H (G, /G,)-1 


n 


1___+ i=n,s,h. 
2(1+v,) 





The modulus of elasticity for the three-phase 
homogeneous fiber-reinforced composite is calculated 
according to Vanin's formula [1]: 


8Gy d- y, y, -V) 


E,, =y E, +0- y Etan  _ 
" 2-y, +X, +-w(%, -DG/G,) 
Sya 1b) 25W Xy, +-y, (x, -DG/G,) 
? |En 8G), xd-y.) +04 y,xG/G,) 
X+y,+d-y,\(G/G,) 
g, = Gite td-wXG/G,) , 
l-y, +(+y,)(G/G,) 
a eS 6 a 
(l—w, )x + (1+ xy, )(G/G, ) 
o l= wR + + y.xG/G,) 
y =| 22E» , Em) We +X +-w."G/G,) 


2") B, 8G _ 2(K-Dd-y,)+(, -I(K-1+4 2y, (G/G,) 
2+xy,—y, +1- y, (x, -1)(G/G, ) 
I+D -v V, 


Va =V- = = 5 
2—w, +xw, +(l-y, (x,-1)(G/G, ) 


E = _ 
Vo =—v,,3xk =3-4V;3 
By 


IV. DOMINANT EQUATION 
4.1. The relationship between strain and displacement 


When taking into account the deformation of the mean 
surface of the element, the strain vector components are 
related to the displacement field according to the 


expression: 
o N 
oy oy K, Es 
= o N 
E, p =3 E, p +Z3 K, 7 +458 ĵo (3) 
o N 
ic Ta Kyy ee 
Where: 
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0 
a Ee 
e Ox R, 
Ex üo 
ji 0 1 . l ; 
e, ¿=| 0 E Vo ¢ is the linear strain vector, 
Yzy 3 ; / Wo 
—— 0 
Oy OX 


Rx, Ry are the radius of curvature in the x and y 
directions, respectively 


T 
{x}= {x K, K = a ain Ogee is the 
n = = Ox Oy Oy Ox 


bending and torsion curvature vectors, 








OW 9 0 
es Ox Oo 
e =| 0 ANa Wo isthe nonlinear strain vector. 
N S: 
Yy ôw, Ow, | LOY 
ôy ê 


4.2. Relationship between stress and strain 
O; Q;; Q,» Qi fy 
Cy ¢ =| Qr Qz Qx E (4) 
Oxy J, Qis Qe Qes Low 

Where: 

Q,, =Q,, cos* 0 + Q,, sinf 6+ 2(Q,,+2Q,.)sin® Ocos*O 

Q,, =Q,,(cos* 0+sin* 0) +(Q,,+Q,,-4Q,,) sin* 8cos*O 


Qi = (Q2 —Q +2 Qs) cos O sin’ 0) +(Q),- Q).- 2 Q,,) in Ocos’0 


Q,, = Q, sinf 6+ Q,,cos*8 + 2(Q,,+2Q,,.)sin* Ocos*O 


Qs = (Qi = em =2 om )cos 0 sin? 0) + (Q, = Q +2 On )sin Q@cos°@ 
Qs = Qes (cos 0 t siní 0) t Qi + Qs _ 2(Q,5 oF Qis )) sin? Ocos70 


E 
= ~ 2 O 
1- (E/E w Yi 
Q, = VQ» Qes = Giz 


4.3. Internal force components 


The surface force vector {N} = {Nx Ny Nxy}', bending 
moment, torsion moment, {M} = {Mx My Myy}! in the 


Ey 2G 7 
11? 


Q, a o Saa 
1-(E, /E Vp Ey, 


shell element with n composite layers are determined as 
follows: 


BHARA  * 


Where: 
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k+ 


mj-( 1o] (a) 


4.4. Differential equation of vibration of a three phase 


rm 
y 
a | 
II 
AN 
n| = 
M: 
~O 
L | 
= 
——~ 
N 
N 
N 
FN 
No 
S v 


composite shell 


After building the finite element model of the structure 
is built, the structural analysis is performed. This work 
includes: 


- Build element equations (element stiffness matrix, 
element load vector); 

- Connect elements to create the overall stiffness 
matrix; 

- Set up the general equation; 

- Solve the general equation; 

- Calculate the necessary results from the solutions of 


the general equation. 


The nonlinear differential equation which describes the 
vibration of the shell: 


[M, ]{a}+[c, ]{a}+[K, ]{a} = {F} (6) 
Where: | M, | - Matrix of the overall mass of the shell; 


| K a - the the overall stiffness matrix of the shell; LC, | 3 


The overall resistance matrix of the shell, calculated by the 
formula: C, | = (a[ m, |+B[ K, }) 


Equation (6) is a nonlinear differential equation that is 
solved by the Newmark method combined with the 
Newton-Raphson iterative method. 


V. NUMERICAL RESULTS AND DISCUSSION 
5.1. The starting problem: 


Structural parameters: Three - phase composite pillar 
panel, rectangular projection size, total thickness h = 
0,0025m, radius of curvature R = 1,0m, length L = 0,30 m, 
opening angle 0 = 30°. The composite shell consists of 5 
layers, the composite layers are made of Graphite/Epoxy 
T300/976, each layer has a thickness of hı = 0,0005m; the 
ratio of grain and matrix is 0,3. Considering the case of 
composite layers arranged symmetrically [-a/ a /0/ a /- al, 
with a = 45°, 

Material characteristics of Graphite/Epoxy T300/976 
are as follows: Graphite-Epoxy T300/976 : Ej; = 150.10? 
N/m’, Ex = E33 = 9.10° N/m’, Gi2 = G13 = 7,1.10° N/m’, 
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Gz = 2,5.10° N/m’, viz = v23 = V32 = 0.3, pce = 1600 the nonlinear method is more stable and safer. According 
kg/m’. to the authors, this is the advantage of solving nonlinear 


Load parameters: a short-term load in the form of a paon ene sou S par t case: 





























shock wave distributed on the upper surface of the shell, a E ee ee eels za 
the load is as follows: | Linear 
p(t) = Pmax FCt), z 4 
= 
t 2 
[==> USTST, (7) - 
F(t) = i . 
0: t >To | 
In which: pmax = 1.105 N/m?, te= 0,025s. The shell was a eae 


sampel along PN Oia enc UE SS = ne ee Fig. 6: Time history response of vertical displacement W 


0, 8, = 0, 8, = 0 at x = O and x = a. The center point on the at point A 


upper surface of the shell (point A (Fig.3)) was considered. 
The finite element model of the problem is shown in Fig.4. 
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Fig. 7:. Time history response of stress ox at point A 


Fig. 3: Real model of the problem 7 
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Fig. 8: Time history response of stress o, at point A 


Fig. 4: Finite element model of the problem 


Table 1. Maximum displacement and stress at point A 


5.2. The influence of some factors on the nonlinear 
vibration of the shelll 


(N/m?) (N/m?) 
5.2.1. Effect of nonlinear properties 


| a | 0,00096 | 4,302.107 | 0,8305.107 
To examine the effect of nonlinearity, the linear problem Cases : - - 
was compared with the solved nonlinear problem. Fig. 6 0,00077 | 3,540.10 0,6470.10 
and Fig.6, Fig.7, Fig.8 and table 1 show the displacement Compare (%) 24,68 21,53 28,36 


and stress variations at point A in two cases. The time 5.2.2. Effects of particles ratio 





response of displacement, stress at the calculation point of 


the linear problem is different from the nonlinear problem To consider the influence of particles ratio, the author 


both in amplitude and cycle. In particular, the response proceeds to solve the problem with 3 cases of different 


values of the nonlinear problem are much larger than that particles ratio. 


of the linear problem, which shows that the calculation by 
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Do vong Wim] 
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Thoi gian t[s] 


Fig. 9: Vertical displacement variation at point A with 
different particles ratio 


The time response of displacement at the shell point in 
the cases of different particles ratio is different in both 
amplitude and period. In which, the deflection response 
values in the case of 0.2 is the largest, in the case of 0.4 is 
the smallest in the 3 cases. This shows that in the 
reasonable range of the ratio between the particles and the 
matrix, when increasing the ratio of the particles, the 
response to the maximum displacement will decrease. 


5.2.3. Effect of load amplitude 


Investigate the problem with load amplitude pmax 
varying from 0.5,105N/m? to 1,66,105N/m°. The results of 
the variation of the maximum values of displacement and 
stress at point A of the shell are shown in graphs in Fig.10, 
Fig.11, Fig.12. 


=] 


#10 


Chuyén vi lön nhat [m] 





%5 LTS 1 235 15 16 
Biên dé täi trong Ppa dem] z 10" 


Fig. 10: Displacement Wmax by amplitude pmax of load 


Ung suat lón nhat Xicma [N/m?] 





$5 0.75 


1 125 15 186 
Biên dô tai trong Pra N/cm?) x 10° 


Fig. 11: Maximum stress by amplitude pmax of load 
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Chuyén vi thang difng tai diém A [m] 


0 apes re as 0.08 0.1 

Thai gian t[s] 
Fig. 12: Time history response of vertical displacement W 
at point A 


The load amplitude has a great influence on the 
dynamic response of the shell. For the survey problem, 
when the load amplitude increases, both the maximum 
displacement and the maximum stress at point A increase 
nonlinearly. The rate of increase of the maximum values of 
displacement and stress is large when pmax is greater than 
1.25.105N/m2. When the load amplitude pmax = 
1,66,105N/m2, both displacement and stress increase very 
strongly, the shell oscillates and balances to another 
position, with hysteresis. According to the dynamic 
stability criterion of Budiansky, B and Roth, R.S [1], the 
shell is destabilized and the critical amplitude value. 


VI. CONCLUSIONS 


This paper presents a finite element method to analyze 
the dynamic response of a 3-phase composite cylindrical 
shell under the action of dynamic loads. The formulations 
are based on classical multilayer shell theory taking into 
account geometric nonlinearity. With the numerical survey 
results on the problem classes with the change of many 
parameters, it is the basis for the selection of reasonable 
parameters for the 3-phase composite cylindrical shell 
structure subjected to dynamic loads. 
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